Finite-size scaling of pseudo-critical point distributions 
in the random transverse-field Ising chain 
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We study the distribution of finite size pseudo-critical points in a one-dimensional random quan- 
tum magnet with a quantum phase transition described by an infinite randomness fixed point. 
Pseudo-critical points are defined in three different ways: the position of the maximum of the av- 
erage entanglement entropy, the scaling behavior of the surface magnetization, and the energy of a 
soft mode. All three lead to a log-normal distribution of the pseudo-critical transverse fields, where 
the width scales as L _1// " with v = 2 and the shift of the average value scales as L~ 1 ^ vtyp with 
vtyp = 1, which we related to the scaling of average and typical quantities in the critical region. 

PACS numbers: 05.50.+q, 75.10.Nr, 75.40.Mg, 75.50.Lk 
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Quenched disorder has a profound effect on the phys- 
ical characteristics of phase transitions in classical and 
quantum mechanical systems. A theoretically and exper- 
imentally important issue is the measurement of physi- 
cal observables in disordered systems at or near critical 
points. These measurements are always performed on 
finite samples and on one particular realization of disor- 
der. Finite size scaling (FSS) pj is the systematic way 
to extract informations on the thermodynamic limit by 
studying finite systems, and the objects to be analyzed by 
FSS of disordered systems are the distributions of physi- 
cal properties in the ensemble of the disorder realizations. 
This also sheds light on the question about whether a sin- 
gle experimental measurement on a rather large system 
is representative for the whole ensemble of random sys- 
tems to which it belongs. This is very much connected 
to the important issue of self-averaging [3, 0] of thermo- 
dynamic quantities like the expectation values for order 
parameter, specific heat or susceptibilities. 

In an infinite system these observables display a char- 
acteristic singularity at a critical point, where for in- 
stance the susceptibility diverges. This divergence is sup- 
pressed in a finite system and replaced by a finite maxi- 
mum, the location of which is called pseudo-critical point 
and is slightly shifted against the infinite system's critical 
point. In pure systems this shift depends in a systematic 
way on the lateral systems size L, usually proportional to 
L~ 1 l Vp , where vp is the correlation length exponent of 
the pure system. In finite disordered systems the suscep- 
tibility usually has several maxima in the critical region 
and each one is slightly shifted against the critical point 
of an infinite system. One identifies the location of the 
largest maximum with the pseudo-critical point of the 
corresponding sample and an intriguing question there- 



fore concerns the distributions of these pseudo-critical 
points. 

If the disorder is irrelevant according to the Harris cri- 
terion vp > d/2 |4|], d being the dimension of the varia- 
tion of the disorder (usually identical with the system's 
spatial dimension), the width of the distribution scales 
as L~ d / 2 . In this case the shift of the average finite-size 
transition point is proportional to L~ 1 l Vp as the shift of 
the finite size pseudo-critical point in the pure case. 

For relevant disorder vp < d/2, there is a new random 
(R) fixed point at which the exponent, vp ^ vp, satisfies 
the relation [5| vp > d/2. In many random systems (di- 
luted magnets, random field problems, spin glasses, etc.) 
the random fixed point is of conventional form, which 
means that thermal and disorder fluctuations remain of 
the same order at large scales, i.e. during renormaliza- 
tion. According to the finite-size scaling theory of con- 
ventional random critical points 3] both the shift and 
the width of the distribution of pseudo-critical points are 
characterized by the same random exponent and there 
is a lack of self-averaging 13]. These predictions were 
checked for various models @, H, [E S @] • 

In this paper we intend to examine the distribution 
of pseudo-critical points and its finite size scaling for 
a quantum phase transition in a paradigmatic model 
for a random magnet with an infinite randomness fixed 
point. This new type of random fixed points has been 
observed in various systems (disordered quantum mag- 
nets at T = 0, random stochastic models, etc.) in which 
the disorder plays a dominant role over quantum, ther- 
mal, or stochastic fluctuations : during renormalization 
the strength of disorder grows without limits [9(. Many 
asymptotically exact results have been obtained for one 
dimensional systems, partially by the use of a strong dis- 
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order renormalization group (SDRG) method [Tjj. 

We study in detail the random transverse-field Ising 
spin chain[ll| (RTFIC) defined by the Hamiltonian: 

ff = -X>ofof +1 -$>af , (1) 

l i 

in terms of the af' z Pauli matrices at site I. Here the Ji 
exchange couplings and the hi transverse-fields are inde- 
pendent random variables. We are interested in the prop- 
erties of the system in its ground state, i.e. at T = 0. In 
the thermodynamic limit the control-parameter is given 

by[H 

5=[ln/i] av -[lnJ] av , (2) 

where [. . . ] av denotes averaging over quenched disorder, 
so that 5 = at the critical point [T^. In the following, 
we consider the case of random coup lings Ji and homo- 
geneous transverse fields hi — h [13J, so h is the analog 
of the temperature in thermal transitions. Its critical 
value in the thermodynamic limit is given by ln/i c (oo) = 
[In J] a y In the vicinity of the critical point the average 
correlation length involves the exponent v = 2, whereas 
the typical correlation length diverges with a different 
exponent v tyv — 1 [9j]. The characteristic time-scale, r, 
which is related to the smallest gap as r ~ e -1 , scales log- 
arithmically at the critical point, lnr ~ y/L. It remains 
divergent also in an extended region of the off-critical re- 
gion, in the so called Griffiths phase, where r ~ L z with 
a dynamical exponent, z, which depends on the distance 
from the critical point. 

Finite-size scaling of the RTFIC has been studied in 
[Til . [TH HH ] . The distribution of the surface magnetiza- 
tion could be computed analytically [HI, [l(| , and the dis- 
tribution of the the gap and the end-to-end correlation 
function by the SDRG method [l6|. They turned out 
to be different in the micro-canonical and the canonical 
ensembles. Here we define sample dependent critical pa- 
rameters h c (a, L) (where a indicates a particular disorder 
realization) and study their distribution. The standard 
approach, defining h c (a,L) through the rounding of the 
singularity of the susceptibility, does not work here, since 
the susceptibility is divergent also in the Griffiths phase. 
A similar problem arises for the rounding of the specific 
heat, since it has only a very weak essential singularity. 

a. Pseudo- critical points through the maximum of the 
average entropy Here we suggest that for random quan- 
tum systems the pseudo-critical transition points can be 
conveniently defined through the rounding of the average 
entanglement entropy. For the RTFIC we consider a pe- 
riodic sample (a) of length L and calculate the entangle- 
ment entropy between the two halves of the chain, which 
is then averaged over all possible starting points of the 
block. In the limit L — > oo the average entropy is diver- 
gent at the critical point (l7j . whereas in a finite sample 
we use the position of its maximum to define h c (a,L). 

In the numerical calculations we have used efficient 
free-fermion techniques [HI, [r| with which we could cal- 
culated the entanglement entropy up to sizes L = 512. 
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FIG. 1: (Color online) Distribution of lnh c (a, L) for a disor- 
der with AJ = 0.4 for different sizes. In the inset the distri- 
butions of scaled variables is well described by a Gaussian. 

The couplings in Eq.([T]) are taken from a uniform box- 
like distribution, which is centered at J = 1 and has a 
width AJ. For different strengths of disorder (AJ = 
0.2, 0.4, 0.6 and 1.0) and for each system size L 10000 
disorder realizations were generated. Additionally the 
entropy of each sample is averaged over L/2 starting po- 
sition of the block. As an illustration we show in FigQ] 
the probability distributions of \nh c {a, L) at a disorder 
AJ = 0.4 for different sizes. The distribution functions 
for different L are symmetric and in terms of rescaled 
variables: (In h c (a, L) — [In h c (a, L)] av )/A In h c (a, L) the 
transformed distributions are well fitted by the same 
Gaussian form, as shown in the inset of FigQ] For dif- 
ferent strength of disorder we have analyzed the shift of 
the average value, [ln/i c (a,L)] av , as well as the standard 
deviation, A In h c (a, L), which are shown in FigfS] In- 
terestingly, the average transition point for a given L is 
practically independent of the strength of disorder and 
corresponds to the value in the pure system. 

Our numerical data are compatible with a FSS form 
for the shift that is given by h c (oo) — h c (L) ~ L~ 2 lnL. 
The scaling of the width of the distributions is found to 
follow Ah c (L) ~ i -1 /", where the exponent is given by 
1/v = 0.50(1) independently of the strength of disorder 
(FigfTJ). Thus, our numerical estimate for v agrees well 
with the exponent of the average correlation length of 
the RTFIC. We can thus conclude that the distribution 
of In h c (a,L) in Fig Q] is well described by a Gaussian 
with a variance of 0(1/ L) and by a shift of 0(L~ 2 lnL). 

b. Pseudo- critical points through the surface magne- 
tization The surface magnetization is perhaps the sim- 
plest physical quantity of the RTFIC. Fixing the spin at 
one end of the chain say at I = L + 1, which amounts 
to have = 0, the magnetization at the other end of 
the chain at / = 1 is given by the exact formula [20j: 

m s = [l + Ef =i ni =1 (^A') 2 ]- 1/2 - (3) 
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FIG. 2: (Color online) Upper panel: Shift of the average of 
In h c (a,L) as a function of L in a log-log scale. The data 
points have practically no disorder dependence. The straight 
line has a slope of 1.8, which represents an effective exponent 
close to 2. Lower panel: Scaling of the width of the distribu- 
tion of In h c (a, L) as a function of L in a log-log scale. The 
slope of the straight lines is practically disorder independent: 
1/^ = 0.50(1) 



In the following we use a doubling procedure |8|: For a 
given random sample (a) of length L, we construct a 
replicated sample (2a) of length 2L by gluing two copies 
of (a) together, and study the ratio of the surface mag- 
netizations: r(a,L) — m s (2a,2L)/m s (a, L). We rewrite 
the exact expression in Eq.© as 



,(a,L) = [Zi(L) + l] 



-1/2 



Z l( L)=X>-^) 
1=1 



(4) 



in terms of a random walk variable: U(i) — 2 Y]l—i In jfr- 
Similarly, we obtain for the replicated sample 

2L 

m s (2a, 2L) = [Z 2 (L) + \y 1/2 , Z 2 (L) = £ e~ u ^ 

i=i 

(5) 

where 1/2(1) = U(i) for 1 < i < L and it is U 2 (i) = 
U(L) + U(i - L) for L + 1 < i < 2L. The expression in 
Eq.© simplifies into Z 2 (L) = (l + e~ u ^) Z X (L). Since 
in the critical region In Z\(V) ~ L 1 / 2 for large L the ratio 
of the two surface magnetizations is given by: 



r(a, L) 



1 + e 



-U(L) 



-1/2 



(0) 



For the pure chain in the ordered phase with Sp — 
h — h c (oo) = h — 1 < we have U(L) — > 00, and in a large 
finite system the ratio behaves as r(L) = 1 — e~ 2 \ Sp \ L /4, 
whereas in the disordered phase with 5p > we have 
U(L) — * — 00 and the finite-size correction reads as 
r(L) = e~ 5pL . Then at the critical point U(L) = 
and the ratio has a non-trivial value, r(L) = l/y/2. For 
the random chain we have the same type of trivial fixed 
points corresponding to the disordered (U(L) — ► 00) and 
the ordered (U(L) — > -co) phases, respectively, and it is 
natural to use the condition, U (L) — 0, to define finite- 
size critical transverse fields. This leads to the micro- 
canonical condition: 



1 L 1 L 

— In Ji = — ^2 m h — ln h c (a, L) 



(7) 



i=i 



Using this definition we obtain for the finite-size be- 
havior in the disordered phase for a given sample 
lnm s (a,L) ~ — (In h c (a,L) — \nh)L and for its average: 
[ln m s (a, L)] av ~ — (ln/i c (oo) — hih)L ~ — \6\L, which 
involves the typical exponent, v typ = 1. 

c. Pseudo-critical points through a soft mode Criti- 
cality of the RTFIC in the free fermion representation is 
related to the vanishing of the excitation energy of a spe- 
cial fermionic mode. We recall that for closed chains, i.e. 
with Jl 7^ the even and odd number of excitations are 
taken from two different sectors and the ground state is 
the vacuum of the even sector, whereas the first excited 
state is in the odd sector and contains one fermion with 
energy, A\. For the pure system its energy is given by 
Af = h — J = Sp, thus it changes sign at the critical 
point. Here we use the same condition, Ai = 0, to define 
criticality a in finite random system, too. This leads to 
the matrix equation [18]: (A + B)*i = with 



Mi 



(A + B) = 



Ji 

h 2 J2 
hs J 3 



flL-l Jl-1 



(8) 



the solution of which is just the micro-canonical condition 
in Eq.(U|). 

In the following we study the scaling behavior of the 
smallest gap e(a, L) in the ordered phase. For this it is 
more convenient to use open chains with L bonds where 
we have the asymptotic expression [T3|: 



e(a, L) ~ m s (a, L)m s (a, L) -^1 



t = l 



(9) 



provided the scaled gap, e(a, L)L, goes to zero. Here 
we take /il+i = hi. m s (a,L) and m s (a, L) denote the 
surface magnetization at the two ends of the chain, which 
are both 0(<5 1//2 ) in the ordered phase. Consequently for 
a given sample, lne(a, L) ~ — (ln/i — ln/i c (a, L))L and 
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for its average: [lne(a, L)] av ~ — (ln/i — hxh c (oo))L ~ 
— <5L. Thus the finite-size correction involves the typical 
exponent, v tyv = 1. 

<i Log-normal distribution We have used three dif- 
ferent methods to determine pseudo-critical points, all 
of which gave coherent results for the pseudo-critical 
point distributions. We found that the distribution of 
x = In h c (a, L) is Gaussian around x = [In h c (a, _L)] av 



P L ( x ^lnh c (a,L)) = ] J^e-^-^ (10) 

Thus A In h c (a, L) = a/ \f~L and the fluctuations are gov- 
erned by the exponent v = 2. The shift of the av- 
erage, ln/i c (oo) — x, is of 0(L~ 2 In L) from the aver- 
age entropy and zero by the other two methods. The 
average of the critical transverse fields is given by: 
[K(a, L)] av = <?e^l 2L = ft c (oo)(l + 2a 2 /L + ...), thus 
[h c (a, L)] m — h c (oo) = \(L) ~ which is consistent 
with v typ — 1. We can thus conclude that Ah c (L) j X(L) 
tends to infinity, which can be taken as a definition of an 
infinite randomness fixed point. 

We now use these results to explain the role of the 
averaged correlation length 5~ 2 . In the disordered 
phase (h > h c (po)), the average surface magnetization 
[m s (h, L)] av is dominated by the rare ordered samples 
having m s (a,L) = O(l), for which h c (a,L) > h > 
h c (oo). From their probability we obtain: 

[m s (h,L)] av ~ Prob(lnh c (a,L) > In h) ~ e~^ 8 ' (11) 



Its exponential decay thus involves the critical exponent 
v = 2. In the ordered phase [h < h c (oo)) we consider 
the average gap, [e(h, £)] a v, which is dominated by those 
rare realizations, for which h c (a,L) < h < h c (oo) and 
e(a, L) = 0(1). From their probability we obtain: 

[e(h, £)] av ~ Prob{\n h c (a, L) < In h) ~ e"^ 5 ' (12) 

which also involves the critical exponent v = 2. We note 
that previously we have shown that the typical quanti- 
ties, both the surface magnetization and the gap, have 
an exponential decay: e _i ' (5 ', which involve the typical 
exponent, Vt yp = 1. 

Our results for the RTFIC are also relevant for other 
systems displaying an infinite randomness fixed point, 
like the random antiferromagnetic XX chain and vari- 
ous stochastic models with quenched disorder, such as 
the Sinai-walk and the partially asymmetric exclusion 
process. Also in higher dimensional realizations of in- 
finite randomness fixed points, like the two-dimensional 
random transverse Ising model [2lL I22I ] , we expect a sim- 
ilar scenario as we described here, which can be checked 
for instance by calculating the entropy numerically by 
the SDRG method H- 
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